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ABSTRACT 

We report the discovery of three extrasolar planets that transit their moderately bright 
(my = 12-13) host stars. WASP-44b is a 0.89-M Jup planet in a 2.42-day orbit around 
a G8V star. WASP-45b is a 1.03-Afj up planet which passes in front of the limb of 
its K2V host star every 3.13 days. Weak Ca II H+K emission seen in the spectra 
of WASP-45 suggests the star is chromospherically active. WASP-46b is a 2.10-Mj up 
planet in a 1.43-day orbit around a G6V star. Rotational modulation of the light 
curves of WASP-46 and weak Ca n H+K emission in its spectra show the star to be 
photospherically and chromospherically active. 

We imposed circular orbits in our analyses as the radial velocity data are consistent 
with (near-) circular orbits, as could be expected from both empirical and tidal-theory 
perspectives for such short-period, ^Jupiter-mass planets. We discuss the impact of 
fitting for eccentric orbits for such planets when not supported by the data. The derived 
planetary and stellar radii depend on the fitted eccentricity and these parameters 
inform intense theoretical efforts concerning tidal circularisation and heating, bulk 
planetary composition and the observed systematic errors in planetary and stellar 
radii. As such, we recommend exercising caution in fitting the orbits of short period, 
^Jupiter-mass planets with an eccentric model when there is no evidence of non- 
circularity. 

Key words: planets and satellites: individual: WASP-44b, WASP-45b, WASP-46b 
- stars: individual: WASP-44, WASP-45, WASP-46 



1 INTRODUCTION 

The ensemble of well-characterised transiting extrasolar 
planets is growing at pace, with well over one hundred 
known to date. It is important to determine the system pa- 
rameters accurately so that the inferrences based on them 
are reliable. For example, to determine the bulk composi- 
tion of a pl anet it is necessary to measure accurately its 



radius (e.g. iFortnev. Mar lev fe Barnes] 2007|). Many short- 
period, giant planets (e .g. WASP- 17b, I Anderson et al.ll20ld . 
I Anderson et al.ll2011bl ) are larger than predicted by stan- 
dard cooling theory of irradiated, gas-giant planets (e.g. 
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IFortnev. Marlev fc Barnesll2007l ). One potential explanation 
is that energy from the tidal circularisation of eccentric or- 
bits was dissipate d within the planets' interiors, ca using 
them to bloat (e.g. iBodenheimer. Lin fc Mardhngll200ll h To 
evaluate the likelihood that a planet was inflated by such 
tidal heating, it is necessary to have an accurate determi- 
nation of both its radius and its orbital eccentricity (e.g. 
Ilbeui. Spiegel fc Burrowsll201lf ). 

A planet's orbital eccentricity can be determined by 
measuri ng the radial motio n of its host star around its or- 
bit fe.g. lQueloz et al.1l2010l h or by observing occultation s of 
the planet by its host star (e.g. I Anderson et a.1.1 l20l"lbh . or 
from a combination of the two. By combining this eccen- 
tricity measurement with high-quality transit light curves, 
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we ca n measure a star's density dSeager fc Mallen-Ornelasl 
2003). The stella r mass can be estimate d using stellar evolu- 
tion mod els (e.g.lDemarque et al.ll2004|) or m ass-calibration 
laws (e.g. iTorres. Andersen fc Gim cncz 2010j), and the stel- 
lar radius follows. This, combined with the ratio of the radii 
derived from the transit depth, gives the planet radius. 

In papers announcing new planet discoveries, eccen- 
tricity is often poorly determined, thus the planet ra- 
dius can be uncertain. D espite this, the re are theoretica l 
l|Goldreich fc Soteri [19661 ) and empirical i|Pont et al.l 120111 ) 
reasons to expect short-period (< 4 d), Jupiter-mass (« 0.5- 
2 Mj U p) planets, often referred to as hot Jupiters, to be in 
circular orbits. Thus it is reasonable, in the absence of evi- 
dence to the contrary, to assume that a newly discovered hot 
Jupiter is in a circular orbit. In so doing, the derived stellar 
and planetary dimensions will be, on the whole, accurate. 

In this paper we present three new hot Jupiters, WASP- 
44b, WASP-45b and WASP-46b, and discuss the effects of 
fitting an eccentric orbit model to a hot Jupiter system in 
the absence of evidence for non-circularity. 
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2 OBSERVATIONS 

2.1 Discovery photometry 

The WASP (Wide Area S earch for Planets) photometric sur- 
vey (|Pollacco et al.ll2006T ) employs two eight-camera arrays, 
each with a field of view of 450 deg 2 , to monitor bright stars 
(my = 8-15) . Each array observes up to eight pointings per 
night with a cadence of 5-10 min, and each pointing is fol- 
lowed fo r around five mont hs at a time. The WASP-South 
station IjHellier et al.l 12011! ) is hosted by the South African 
Astr onomical Observato ry and the SuperWASP-North sta- 
tion (|Faedi et al.ll2011bT ) is hosted by the Observatorio del 
Roque de Los Muchachos on La Palma. The WASP data 
w ere processed and searched fo r transit signals as described 
in lCollier Cameron et ah d2006T) and the candidate s election 
process is described in lcollier Cameron et al.l (|2007t ). 

WASP-44 is a m v = 12.9, G8V star located in the con- 
stellation Cetus. A transit search of the light curve obtained 
by WASP-South from 2009 July to November found a weak, 
2.42-d periodicity. Further observations in 2008 and 2009 
with both WASP instruments led to a total of 15 755 pho- 
tometric measurements (Figure [TJ). WASP-44 was also ob- 
served by a SuperWASP-North camera during 2010 August 
to November. The resulting light curve comprises 6 000 mea- 
surements and became available during the preparation of 
this paper. As such, it was not used in determination of the 
system paramters (Section 3]), but was used in the search for 
rotational modulation (Section [5]). 

WASP-45 is a m v = 12.0, K2V star located in the con- 
sellation Sculptor that was observed by WASP-South during 
May to November of 2006 and 2007. A transit search of the 
resulting light curve, which comprises 11007 photometric 
measurements, found a strong 3. 13-d periodicity (Figure [2]). 

WASP-46 is a m v = 12.9, G6V star located in the con- 
stellation Indus that was observed by two WASP-South cam- 
eras during May to October of 2008 and 2009. A transit 
search of the resulting light curve, which comprises 41 961 
photometric measurements, found a strong 1.43-d periodic- 
ity (Figure [3]). WASP-46 was also observed by two WASP- 



Figure 1. Photometry of WASP-44, with the best- fitting transit 
model superimposed. Top: WASP discovery light curve, folded 
on the ephemeris of Table [6] Middle: WASP data around the 
transit, binned in time with a bin width of ~9 min. Bottom: 
High-precision transit light curve from Euler. 



South cameras during 2010 August to November. The result- 
ing light curve comprises 19 000 measurements and became 
available during the preparation of this paper. As such, it 
was not used in the determination of the system paramters 
(Section [4]), but was used in the search for rotational mod- 
ulation (Section [5}. 

WASP-46 is blended in the WASP images with a m v 
= 14.8 star at a separation of 17.4" (the plate scale is 
13.7"pixel _1 and the photometry aperture has a radius of 
3.5 pixels). We corrected the WASP photometry for this 
contamination to prevent dilution of the transit. 



2.2 Spectroscopic follow-up 

In 2010 we used the CORALIE sp ectrograph mounted 
on the 1.2-m Eule r-Swiss telescope (jBaranne et al.l 1 19961 : 
lOueloz et al. 2000) to obtain spectra of the three target 
stars. We obtained 17 spectra of WASP-44, 13 spectra of 
WASP-45 and 16 spectra of WASP-46. Radial-velocity (RV) 
measurements were c o mputed by w eighted cross-correlation 
l|Baranne et al.l 1 19961 ; iPepe et al.l 120051 ) with a numerical 
G2-spectral template for both WASP-44 and WASP-46, 
and with a K5-spectral template for WASP-45. RV vari- 
ations were detected with the same periods found from 
the WASP photometry and with semi-amplitudes consistent 
with planetary-mass companions. The RV measurements are 
listed in Table [1] and are plotted separately for each system 
in Figures H [5] and H We excluded two RVs of WASP-44 
from the analysis as the spectra were taken during transit 
and we do not fit f or the Rossiter-McLaughlin effect (e.g. 
ISimpson et~al]|201lD . 

For each system we tested the hypothesis that the RV 
variations are due to spectral-line distortions caused by 
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Figure 2. Photometry of WASP-45, with the best-fitting transit 
model superimposed. Top: WASP discovery light curve, folded 
on the ephemeris of Table [6] Middle: WASP data around the 
transit, binned in time with a bin width of ~11 min. Bottom: 
High-precision transit light curves from Euler (upper) and TRAP- 
PIST (lower). 
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Figure 3. Photometry of WASP-46, with the best-fitting transit 
model superimposed. Top: WASP discovery light curve, folded 
on the ephemeris of Table [6] Middle: WASP data around the 
transit, binned in time with a bin width of ~5 min. Bottom: 
High-precision transit light curves from TRAPPIST (upper) and 
Euler (lower). 



a blended eclipsing binary or starspots by performing a 
line-bisector analysis l|Queloz et al.ll200ll ) of the CORALIE 
cross-correlation functions. The lack of correlation between 
bisector span and RV (Figures 01 [S] and [S| supports our con- 
clusion that the periodic dimming and RV variation of each 
system are caused by a transiting planet. 



2.3 Follow-up photometry 

We obtained high signal-to- noise (S/N) transit photometry 
to refine the parameters of each system. We used the 1.2-m 
Euler-Swiss telescope located at ESO La Silla Observatory, 
Chile to observe one transit each of WASP-44b, WASP-45b 
and WASP-46b. Observations were made through a Gunn r 
filter and Euler's absolute tracking mode was used to keep 
the stars on the same location on the chip. This is done by 
calculating an astrometric solution for each science frame 
and adjusting the telescope pointing between exposures. Af- 
ter correcting the images for bias, overscan and flat-field 
variations, we extracted aperture photometry for the bright 
stars in the field. We selected the aperture radius and set 



of reference stars that resulted in the smallest light curve 
residuals. 

We used TRAPPIST U, a 60 cm roboti c telescope also 
located at ESO La Silla (foillon et al.ll201ll ) to observe one 
transit of WASP-45b and one transit of WASP-46b. Both 
transits were observed through a special I+z filter, in the 
2 MHz read-out mode, and with lxl binning, resulting in 
a typical combined readout and overhead time of 8 s and a 
readout noise of 13.5 e~ . The 'software guiding' system of 
TRAPPIST kept the stars at the same positions on the chip 
to within a few pixels over the course of the observations 
(Gillon et al., in preparation). 

2.3.1 WASP-Ub 

A transit of WASP-44b was observed on 2010 September 14 
by Euler. We observed for a total of 4.2 hours, from 46 min 
before the start of the transit to 65 min after it ended. The 

1 TRAnsiting Planets and Planeteslmals Small Telescope; 
http://arachnos.astro.ulg.ac.be/Sci/Trappist 
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Table 1. Radial velocity measurements 



Star 


BJD (UTC) 


RV 


0"RV 


BS 




(day) 


(km s J 


(km s ) 


(km s J 


WASP-44 


2 455 378.8533 


-4.114 


0.028 


-0.063 


WASP-44 


2 455 387.7850 


-3.923 


0.019 


-0.011 


WASP-44 


2 455 388.9089 


-4.173 


0.024 


-0.020 


WASP-44 


2 455 389.7778 


-3.894 


0.040 


-0.191 


WASP-44 


2 455 390.7802 a 


-4.074 


0.025 


0.006 


WASP-44 


2 455 391.8691 


-4.095 


0.020 


-0.156 


WASP-44 


2 455 392.7604 


-3.905 


0.023 


-0.038 


WASP-44 


2 455 396.8995 


-4.033 


0.039 


-0.013 


WASP-44 


2 455 414.8075 


-3.981 


0.020 


-0.026 


WASP-44 


2 455 446.8271 


-4.149 


0.027 


-0.065 


WASP-44 


2 455 448.9128 a 


-4.013 


0.032 


-0.005 


WASP-44 


2 455 449.8568 


-4.173 


0.024 


-0.012 


WASP-44 


2 455 450.8095 


-3.912 


0.025 


0.001 


WASP-44 


2 455 451.8089 


-4.162 


0.025 


-0.059 


WASP-44 


2 455 453.6439 


-3.999 


0.018 


-0.017 


WASP-44 


2 455 454.8482 


—4.049 


0.026 


-0.069 


WASP-44 


2 455 482.6600 


-3.978 


0.025 


-0.096 


WASP-45 


2 455 388.9353 


4.406 


0.010 


-0.044 


WASP-45 


2 455 390.8149 


4.680 


0.012 


0.030 


WASP-45 


2 455 391.8937 


4.414 


0.011 


0.016 


WASP-45 


2 455 392.8314 


4.552 


0.010 


0.012 


WASP-45 


2 455 393.7593 


4.704 


0.014 


0.013 


WASP-45 


2 455 404.8466 


4.426 


0.010 


-0.035 


WASP-45 


2 455 410.7215 


4.392 


0.019 


-0.048 


WASP-45 


2 455 414.7627 


4.557 


0.015 


-0.011 


WASP-45 


2 455 446.8614 


4.683 


0.012 


-0.024 


WASP-45 


2 455 449.8794 


4.699 


0.011 


0.018 


WASP-45 


2 455 451.8550 


4.452 


0.015 


-0.000 


WASP-45 


2 455 453.8741 


4.507 


0.011 


-0.004 


WASP-45 


2 455 482.6826 


4.410 


0.013 


0.033 


WASP-46 


2 455 334.8428 


-3.414 


0.029 


-0.044 


WASP-46 


2 455 359.8870 


-4.116 


0.033 


-0.014 


WASP-46 


2 455 382.7544 


-4.092 


0.062 


0.158 


WASP-46 


2 455 385.8357 


-3.869 


0.032 


-0.016 


WASP-46 


2 455 389.8330 


-4.144 


0.024 


-0.075 


WASF-46 


2 455 390.0889 


—3.491 


0.036 


—0.074 


WASP-46 


2 455 391.7867 


-3.516 


0.026 


0.057 


WASP-46 


2 455 392.6562 


-4.183 


0.031 


0.010 


WASP-46 


2 455 396.7744 


-4.058 


0.026 


-0.033 


WASP-46 


2 455 400.6822 


-3.382 


0.027 


0.009 


WASP-46 


2 455 409.6251 


-4.025 


0.028 


0.005 


WASP-46 


2 455 413.5695 


-3.555 


0.033 


0.082 


WASP-46 


2 455 414.8452 


-3.433 


0.047 


0.115 


WASP-46 


2 455 445.6578 


-4.141 


0.024 


-0.092 


WASP-46 


2 455 449.7606 


-4.090 


0.036 


0.047 


WASP-46 


2 455 454.7010 


-3.513 


0.026 


-0.083 



a These two RVs were excluded 
measured during transit. 



from the analysis as they were 



conditions were clear, airmass ranged from 1.04 to 1.44 and 
no defocus was applied. The WASP-44b transit light curve 
is shown in Figure [1] and the data are given in Table [2] 



2.3.2 WASP-45b 

A partial transit of WASP-45b was observed on 2010 
September 11 by Euler. We observed for a total of 3.3 hours, 
from part- way through the transit to 2 hours after it ended. 
The seeing was poor (1-2") and airmass ranged from 1.01 
to 1.22. We applied a defocus of 0.1 mm to increase the duty 
cycle and to minimise pixel-to-pixel effects. 

A transit of WASP-45b was observed on 2010 December 
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Figure 4. Top: CORALIE RVs of WASP-44 with the best- 
fitting circular Keplerian orbit superimposed (rms=25.0ms , 
X 2 =10.9). The RVs represented by open, red circles were excluded 
from the analysis as they were taken during transit. Middle: Bi- 
sector span with respect to RV. We adopted uncertainties on the 
bisector spans twice the size of those on the RVs. Bottom: Resid- 
uals of the RVs about the fit as a function of time. 



17 by TRAPPIST. Observations lasted from 0hl5 to 3h40 
UT, during which the airmass increased from 1.03 to 1.73 
and the transparency was good. The telescope was strongly 
defocused to average pixel-to-pixel sensitivity variations and 
to increase the duty cycle, resulting in a typical full width 
at half- maximum of the stellar images of ~9 pixels (~5.8"). 
The integration time was 6 s. After a standard pre-reduction 
(bias, dark and flat-field correction), the stellar fluxes were 
extracted from the image s using the I RAF/DAOPHCnff] aper- 
ture photometry software (|Stetsonlll987l ). We tested several 
sets of reduction parameters and chose the set that gave the 
most precise photometry for the stars of similar brightness 
to WASP-45. We carefully selected a set of reference stars 
and then performed differential photometry. 

The WASP-45b transit light curves are shown in Fig- 
ure [2] and the data are given in Table [2] 



2 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 5. Top: CORALIE RVs of WASP-45 with the best-fitting 
circular Kcplcrian orbit superimposed (rms=9.7ms _1 , x 2 =7.4). 
Middle: Bisector span with respect to RV. We adopted uncer- 
tainties on the bisector spans twice the size of those on the RVs. 
Bottom: Residuals of the RVs about the fit as a function of time. 



2.3.3 WASP-46b 

A transit of WASP-46b was observed on 2010 July 19 by 
TRAPPIST. Observations lasted from lhl5 to 4h05 UT, 
during which the airmass decreased from 1.83 to 1.22 and the 
transparency was good. Considering the relative faintness of 
the target, we chose to keep the telescope focused during 
the run, and we used an integration time of 60 s. The pre- 
reduction and reduction procedures were similar to those 
used for the WASP-45 data. 

A transit of WASP-46b was observed on 2010 Septem- 
ber 10 by Euler. We observed for a total of 3.4 hours, from 35 
min before the start of transit to 75 min after it ended. The 
conditions were variable and airmass increased from 1.12 to 
1.4 during the observations. No defocus was applied. 

The WASP-46b transit light curves are shown in Fig- 
ure [3] and the data are given in Table [2] 



3 STELLAR PARAMETERS FROM SPECTRA 

For each system, the CORALIE spectra were co-added to 
produce a single spectrum with a typical S/N of around 
50:1. The analysis was performed using the methods given 
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Figure 6. Top: CORALIE RVs of WASP-46 with the best- 
fitting circular Keplerian orbit superimposed (rms=34.2ms , 
X 2 =19.2). Middle: Bisector span with respect to RV. Wc adopted 
uncertainties on the bisector spans twice the size of those on the 
RVs. Bottom: Residuals of the RVs about the fit as a function 
of time. 



Table 2. Euler and TRAPPIST photometry of the three stars 



Set Star 



Tel. 



BJD(UTC) 
-2 450 000 



Rel.flux ARel.flux 



1 


WASP-44 


Euler 


5453.686865 


0.9989 


0.0015 


1 


WASP-44 


Euler 


5453.687942 


1.0007 


0.0015 


1 


WASP-44 


Euler 


5453.689683 


1.0000 


0.0015 


5 


WASP-46 


Euler 


5449.614796 


0.9991 


0.0014 


5 


WASP-46 


Euler 


5449.616092 


0.9978 


0.0014 


5 


WASP-46 


Euler 


5449.617396 


1.0003 


0.0014 



This table is available in its entirety in the online journal. 



in lGillon et al.l (|2009h . The Ha line was used to determine 
the effective temperature (T e s), while the Na I D and Mg I 
b lines were used as surface gravity (log p.) diagnostics. The 
parameters obtained from the analysis are given in Table [3] 
The elemental abundances were determined from equivalent- 
width measurements of several clean and unblended lines. 
A value for microtur bulence (£ t ) w as determined from Fe I 
using the method of iMagaml |l984f) . The quoted error esti- 
mates include that given by the uncertainties in T e e, logg. 
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Table 3. Stellar parameters from spectra 



Parameter (Unit) 


WASP-44 


WASP-45 


WASP-46 


T e3 (K) 


5400 ± 150 


5100 ± 200 


5600 ± 150 


log q* fcesl 


4.5 ± 0.2 


4.4 ±0.2 


4.4 ±0.2 


ft fkms" 1 ) 


1.0 ± 0.2 


0.5 ± 0.3 


1.0 ± 0.2 


v sin I (kms - ^) 


3.2 ± 0.9 


2.3 ± 0.7 


1.9 ± 1.2 


[Fe /HI 

L / J 


0.06 ± 0.10 


0.36 ± 0.12 


-0.37 ± 0.13 


[Na/H] 


0.24 ± 0.09 


0.56 ± 0.14 


-0.35 ± 0.09 


Mg/Hl 

L &/ J 


0.25 ± 0.07 


0.57 ±0.09 


-0.10 ± 0.10 


Al/Hl 

L / J 


0.20 ± 0.13 


0.45 ± 0.11 


-0.24 ± 0.08 


Si/Hi 


0.25 ± 0.13 


0.50 ± 0.15 


-0.27 ± 0.10 


[Ca/H] 


0.10 ± 0.17 


0.38 ± 0.17 


-0.21 ± 0.13 


[Sc/H] 


0.30 ±0.17 


0.45 ±0.27 


-0.13 ±0.16 


[Ti/H] 


0.12 ±0.10 


0.46 ±0.20 


-0.29 ±0.18 


[V/HJ 


n o a _i_ n 1 1 
U.z4 ± U.J.7 


U.Yo ± U.lo 


n o 1 _i_ c\ o 1 


/ TTi 
[Cr/HJ 


U.lo ± U.U/ 


U.v54 ± 0.12 


— O.ol ± 0.12 


[Mn/ri] 


n on I n 1 n 

0.29 ± 0.10 


0.93 ± 0.24 


—0.47 ± 0.14 


[Co/HJ 


0.23 ± 0.09 


U.7U ± U.ll 


— 0.32 ± 0.24 


[Ni/HJ 


0.13 ± 0.12 


O.ol ± O.lo 


n qo _l_ n in 
— 0.39 ± 0.10 


log _/HL;l J 


< 1.0 


< 0.8 


< 0.8 


M„ (Mo) 


0.95 ± 0.08 


0.95 ± 0.10 


0.93 ± 0.09 


R* (-R©) 


0.90 ±0.22 


1.00 ±0.25 


1.00 ±0.26 


Spec. Type 


G8V 


K2V 


G6V 


Gyro. Age 


nn+i.o 

u - M -0.6 


I 4+2-0 


1 4+ - 4 


R.A. (J2000) 


00 h 15 m 36?76 


00 h 20 m 56!99 


21 h 14 m 56!86 


Dec. (J2000) 


-ll°56'17'/4 


-35°59'53'.'8 


-55° 52' 18'.' 1 


m v 


12.9 


12.0 ±0.2 


12.9 


m., b 


11.70 ±0.02 


10.75 ± 0.02 


11.75 ± 0.02 


m H b 


11.41 ±0.03 


10.37 ±0.03 


11.47 ±0.03 


m K b 


11.34 ±0.03 


10.29 ±0.02 


11.40 ±0.03 


2MASS b 


00153675-1156172 


00205699-3559537 


21145687-5552184 



a log A(Li) = logA r (Li/H) + 12, where Af( Li/H]) is the number d ensity of lithium with respect to hydrogen. 

b [Skrutskie et al.l i2006r i. 



and <J t , as well as the scatter due to measurement and atomic 
data uncertainties. 

The projected stellar rotation velocities (v sin I) were 
determined by fitting the profiles of several unblended Fe I 
lines. For this we used an instrumental FWHM of 0.11 ± 0.01 
A (= 5.2 ± 0.5 kms -1 ), determined from the telluric lines 
around 6300 A (equating to a spec tral resolution of ~57 000) , 
and used the lBruntt et al.l (|2010l ) calibration to assume val- 
ues for macroturbulence (w ma c). The « mac values we assumed 
were 1.4 ± 0.3 kms" 1 for WASP-44, 0.7 ± 0.3 kms" 1 for 
WASP-45, and 2.0 ± 0.3 kms" 1 for WASP-46. 

Though the S/N are very low in the Ca II H+K regions 
of our spectra, we searched for signs of emission, which would 
be indicative of stellar chromospheric activity. The WASP- 
44 spectra show no signs of emission and both the WASP-45 
and WASP-46 spectra show weak emission. 

We input our values of T e g, log.q* and [ F e/H1 into the 
calibrations of iTorres. Andersen fe Gimene j (|2010T ) to ob- 
tain estimates of the stellar mass and radius. 

We used the method of I Barnes] |2007t) to estimate gy- 
rochronological ages for the stars. For WASP-44 and WASP- 
45 we calculated rotation periods from our v sin I determina- 
tions. If either star is inclined with respect to the sky plane 
then it will be rotating faster than suggested by vsinl and 
so is likely to be younger than indicated. For WASP-46, the 
rotation period found from the modulation of the WASP 



light curves was used (see Section As such, the age esti- 
mate for this star is not just an upper limit. 



4 SYSTEM PARAMETERS FROM RV AND 
TRANSIT DATA 

4.1 WASP-44b and WASP-46b 

We determined the parameters of each system from a si- 
multaneous fit to all data. The fit was performed us- 
ing the current version of t he Markov-chain Monte Carlo 
(M CMC) code described by ICollier Cameron et all (|2007t l 
and iPollacco et all (|2008l ). 

The transit light curves wer e modelled using the formu- 
lation of lMandel fe Agoll (|2002r i with the assumption that 
the planet is much smaller than the star. Limb-darkening 
was accounted for using a four-coefficient, nonlinear limb- 
darkening model, using coeffi cients appropriate t o the pass- 
bands from the tabulations of lClaretl (|2000l , l2004h . The coef- 
ficients were interpolated once using the values of log and 
[Fe/H] in Table El but were interpolated at each MCMC 
step using the latest value of T c g . The coefficient values cor- 
responding to the best-fitting value of T e g are given in Ta- 
ble [4] The transit light curve is parameterised by the epoch 
of mid-transit To, the orbital period P, the planet-to-star 
area ratio (R p \ / R*) 2 , the approximate duration of the tran- 
sit from initial to final contact T14, and the impact param- 
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eter b = a cos i/R* (the distance, in fractional stellar radii, 
of the transit chord from the star's centre in the case of a 
circular orbit), where a is the semimajor axis and i is the 
inclination of the orbital plane with respect to the sky plane. 

The eccentric Keplerian radial-velocity orbit is param- 
eterised by the stellar reflex velocity semi-amplitude Ki, 
the systemic velocity 7, and -yecosoj and -yesinaj (Collier 
Cameron, in preparation), where e is orbital eccentricity and 
uj is the argument of periastron. 

The linear scale of the system depends on the orbital 
separation a which, through Kepler's third law, depends on 
the stellar mass M*. At each step in the Markov chain, the 
latest values of p*, T e ff a nd [Fe/H] are input in to the em- 
pirical mass calibration of lEnoch et ail fcOloT l to obtain A/* . 
The shapes of the tr ansit light curves and the radial- velocity 
curve constrain p, i|Seager fe Mallen-Ornelasl |2003( ) ■ which 
combines with M* to give 7?*. T ff and [Fe/H] are proposal 
parameters constrained by Gaussian priors with mean val- 
ues and variances derived directly from the stellar spectra 
(see Section [3]). 

As the planet-star area ratio is determined from the 
measured transit depth, R p \ follows from R t . The planet 
mass M p i is calculated from the the measured value of K\ 
and the value of M*; the planetary density p p \ and surface 
gravity log g p \ then follow. We calculate the planetary equi- 
librium temperature T Gq i, assuming zero albedo and efficient 
redistribution of heat from the planet's presumed permanent 
day side to its night side. We also calculate the durations 
of transit ingress (T12) and egress (T34), the total occulta- 
tion duration (Tgs), and the durations of occultation ingress 
(T 56 ) and egress (T 78 ). 

At each step in the MCMC procedure, model transit 
light curves and radial velocity curves are computed from 
the proposal parameter values, which are perturbed from 
the previous values by a small, random amount. The x 2 
statistic is used to judge the goodness of fit of these models 
to the data and a step is accepted if \ 2 is lower than for 
the previous step. A step with higher \ 2 is accepted with 
a probability proportional to exp(— A\ 2 /2), which gives the 
procedure some robustness against local minima and leads 
to the thorough exploration of the parameter space around 
the best-fitting solution. To give proper weighting to each 
photometry data set, the uncertainties were scaled at the 
start of the MCMC so as to obtain a photometric reduced- 
X 2 of unity. For the same reason, a jitter term of 10.2 ms" 1 
was added in quadrature to the formal errors of the WASP- 
46 RVs. A possible source of this jitter is stellar activity, 
as indicated by chromospheric and photometric indicators 
(Sections [3] and [5]). With reduced spectroscopic \ 2 values of 
less than unity, it was not necessary to add any jitter to the 
RVs of WASP-44. 

4.2 WASP-45b 

The analysis for WASP-45b was the same as for WASP- 
44b and WASP-46b, except we used evolutionary models to 
impose priors on stellar mass and density for reasons to be 
explained. 

The apparent lack of second and third contact points 
in the transit light curves of WASP-45b (Figure [2]) indi- 
cate that the transit is grazing or near-grazing. When b ± 
R p \/R,> 1, the planet does not fully pass in front of the 
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Fi gure 7. Evolutionary diagram for WASP-45 with the models 
of iDemarque et al.l l|2004h . The solid, black lines are isochrones, 
with the ages coming from gyrochronology: [Fe/H] = 0.48, 3.4 
Gyr (top); [Fe/H] = 0.36, 1.4 Gyr (middle); [Fe/H] = 0.24, 0.4 
Gyr (bottom). The two dashed, blue lines are the mass tracks 
for 1.0 Mq with [Fe/H] = 0.48 and 0.85 M Q with [Fe/H] = 0.24. 
The isochrones and mass tracks we chose to display are those that 
bound the region of interest, which is the area shaded dark grey 
that indicates the range of stellar densities consistent with the 
models. The solid, red line indicates our best-fitting T e g, and the 
area shaded light grey corresponds to the 1-ct confidence interval. 

star. Thus, as b tends toward higher values, R p \ must inflate 
so that the ratio R p i/R* remains consistent with the ob- 
served transit depth. Also, R* must inflate so as to remain 
consistent with the observed transit duration. This results in 
stellar and planetary radii that are non-physical in a sizeable 
portion of accepted MCMC steps. 

To avoid biasing the derived stellar and planetary radii 
to larger values, we imposed a prior on stellar density us- 
ing the constraints provided by our gyrochronological age 
determination (Table [3]) and evolutionary models. Our gy- 
rochronolgical age is determined using v sin / and so is an 
upper limit. If the stellar spin axis is inclined with respect 
to the sky plane then the star will truly be rotating faster 
and will theref ore be younger. We plot ted the isochrones and 
mass tracks of lDemarque et all (|2004l ) relevant to the value 
ranges we obtained for T e $, [Fe/H] and age (Figure [7]) . By 
interpolating the models for these measured values, we ob- 
tained stellar dimensions consistent with the models: p, = 
1.65 ±0.40 p and M, = 0.91 ±0.06 M©. We used these val- 
ues to place Gaussian priors, by means of Bayesian penalties 
on x 2 , at each MCMC step. 

It was not necessary to add any jitter to the RVs of 
WASP-45. 

4.3 The imposition of circular orbits 

For each system, the best-fitting eccentricity is small and 
co nsistent with zero (Tabl e [3). We used the F-test approach 
of iLucv fe Sweeney! (|l97ll ) to calculate the probability, Pls, 
that the improvement in the fit that results from fitting for 
an eccentric orbit could have arisen by chance if the under- 
lying orbit were circular. For each system, we found Pls to 
be much higher (Table O than t he threshold of 0.05 sug- 
gested by ILucv fe Sweeney! (|l97ll) , below which they con- 
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Table 4. Limb-darkening coefficients 



Planet 


Instrument 


Observation bands 


Claret band 


a\ 


«2 


0-3 


a 4 


WASP-44 


WASP / Eulcr 


Broad (400-700 nm) / Gunn r 


Cousins R 


0.646 


-0.437 


1.080 


-0.519 


WASP-45 


WASP / Euler 


Broad (400-700 nm) / Gunn r 


Cousins R 


0.755 


-0.768 


1.450 


-0.623 


WASP-45 


TRAPPIST 


Cousins Z+Sloan z' 


Sloan z' 


0.833 


-0.891 


1.304 


-0.547 


WASP-46 


WASP / Euler 


Broad (400-700 nm) / Gunn r 


Cousins R 


0.543 


-0.055 


0.592 


-0.348 


WASP-46 


TRAPPIST 


Cousins /+Sloan z' 


Sloan z' 


0.630 


-0.371 


0.727 


-0.366 



sider the detection of a non-zero eccentricity to be signifi- 
cant. In the absence of conclusive evidence to the contrary, 
we assumed the orbit of each planet to be circular in pro- 
ducing the adopted solutions that we present. This issue is 
dealt with at greater length in the discussion. 

4.4 Results 

The median values and 1-cr uncertainties of the system pa- 
rameters derived from the MCMC model fits are presented 
in Table The corresponding transit models are superim- 
posed on the transit photometry in Figures [1] [2] and [3] The 
corresponding orbit models are superimposed on the radial 
velocities in Figures 0] [S] and [51 in which we also show the 
RV residuals as a function of time. From the residual plots 
one can evaluate both the level to which stellar activity may 
have affected the RVs and the evidence for a third body. The 
RVs of each planet-hosting star, obtained over the course of 
100 days, show no evidence for motion induced by a third 
body. The residuals of WASP-46 show greater scatter (x 2 
= 19.2 for WASP-46 compared to 10.9 for WASP-44 and 
7.4 for WASP-45), which is probably due to the star being 
active (see Sections [3] and 0) . 

We also performed MCMCs fitting an eccentric model 
to the RVs so as to illustrate the impact on those parame- 
ters that are liable to be affected when imposing a circular 
orbit. The best-fitting values and associated uncertainties 
from MCMCs fitting an eccentric model are given in Ta- 
ble [S] and these are to be compared with the values of our 
adopted solutions in Table[6j derived using a circular model. 
For each system, the best-fitting eccentricity is small and the 
improvement in the fit is slight. The parameter values of the 
circular and eccentric models agree to well within 1 a, and 
the associated uncertainties of the eccentric model are only 
fractionally larger. 



5 ROTATIONAL MODULATION 

We analysed the WASP light curves of each star to deter- 
mine whether they show periodic modulation due to the 
combination of magnetic activity and stellar rotation. We 
used t he sine- wave fitting method described in lMaxted et alj 
(2011) to calculate periodograms such as those shown in the 
upper panels of Figure [8] These are calculated over 4 096 
uniformly spaced frequencies from to 1.5 cycles day -1 . 
The false alarm probability levels shown in these figures are 
calculate d using a boot-strap Monte Carlo method also de- 
scribed in M axted et alj (|201lT ). Variability due to star spots 
is not expected to be coherent on long timescales as a con- 
sequence of the finite lifetime of star-spots and differential 



rotation in the photosphere so we analysed each season of 
data separately. 

We found no evidence of rotational modulation in any 
of the WASP-44 or WASP-45 light curves. 

We applied the sine-wave fitting method to the WASP- 
46 light curves from each of the three seasons to produce the 
periodograms in the upper panels of Figure [8] The peaks 
in the periodograms for the first two seasons of data are 
highly significant and imply periods of 16.55 ± 0.10 d and 
14.92 ± 0.07 d. We used a boot-strap Monte Carlo method 
to estimate the errors on these periods. The periodogram 
for the third season of data shows peaks corresponding to 
a period of 8.08 d. Our boot-strap Monte Carlo method 
fails to provide a reliable error estimate for this less signifi- 
cant peak. The amplitudes of the variations estimated from 
the sine-wave fits are in the range of 3-5 milli-magnitudes 
(Fig. El lower panels). We also calculated periodograms of 
two stars with similar magnitudes and colours to WASP- 
46 observed simultaneously with the same camera. Neither 
of these stars showed significant periodic modulation in the 
same frequency interval with the exception of one season of 
data for one star that showed significant power close to 1 
cycle day - , the source of which is likely to be diurnal. 

Our interpretation of these results is that the WASP 
light curve of WASP-46 does show periodic modulation due 
to the rotation of the star with a period of 16 ± 1 days. The 
period of 8.08 days derived from the third season of data 
can be explained by the distribution of star spots during 
this season resulting in a light curve with two minima per 
cycle. Similarly, the poorly defined peak in the periodogram 
for the second season of data may be due to the appearance 
and disappearance of star spot groups during the observing 
season. 

Using the rotational-modulation period of 16 ± 1 days 
and the stellar radius from Table [6j we calculate the stel- 
lar rotation velocity v to be 2.9 ± 0.2 kms -1 . That this 
value is higher than the value of wsin/ (1.9 ± 1.2 kms -1 ) 
determined from spectral line profiles (Section [3j suggests 
that the stellar spin axis is inclined by I = 41^2? ° with 
respect to the sky plane. This should be appreciated when 
assessing whether a measurement of WASP-46b's Rossiter- 
McLaughl in effect is indicat ive of an aligned or a misaligned 
orbit fe.g. lSchlaufmanlfeoiOi ). Usually only the sky-projected 
angle between the planetary orbit and stellar spin axis can 
be determined, but with a measurement of / the true an- 
gle can be determined. Our measurement of I is imprecise 
and would be improved by a more precise determination of 
spectroscopic v sin /. 
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Table 5. The limits placed on orbital eccentricity and the impact of fitting for it 



Parameter (Unit) WASP-44 WASP-45 



WASP-46 



rmsspcc (ms x ) 

Xspcc 

e 

Pls 

e (2 a upper limit) 
e (3 cr upper limit) 

U n 

e cos u> 
es'mui 

0mid — occultation 

T 58 (d) 

T 5 e « T 78 (d) 

M, (M Q ) 
R* (R@) 
P* 0©) 

M p i (M Jup ) 
R P i (Rjup) 

Ppl (P.Tup) 



24.8 
10.5 



0.036 



+0.054 
0.025 
0.51 



< 0.172 

< 0.263 
-59lf 7 2 

0.003±°;°o?l 
-o mi+ ' 028 

u.uii_ 069 

c i02 +0 ' 015 

091 ^+ - 0049 
u.uyio_ 0093 

ni32+ - 0033 
u.uraz_ 0030 

0.948 ± 0.034 

n oqn+0.085 
u.o»^_ og7 

,+0.54 



1.33" 



0.30 



n oqo + 0.071 

U.O»0_g 

-,+0.13 
-0.14 



Log"* 



8.8 
6.3 



0.023 



+0.027 
-0.016 
1 

< 0.095 

< 0.148 

9-1+59 

oio +0 - 015 
u.uiu_ 011 

o no3+ 034 

0.5061 
0.0711 
> 0.0257 



+0.0096 
0.0067 

+0.0055 
0.0073 



0.910 ±0.060 
Q50+ 093 

U.9OU_ g 74 

1.06 ±0.27 



0.69 



+0.37 
-0.20 



1.005 ± 0.053 

-, 17 +0.28 
± - 1 '-0.14 

0.64 ±0.30 



33.9 
18.1 

0.018_ 013 
0.67 
< 0.065 



< 0.092 

50 
122 



o.oort° : °s 

0.008±r ?I 

U - OUU '~0.0066 

n n7nns+ 00089 

U - U ' UUS -0. 00104 

O.U184_ 001g 

0.957 ±0.034 
0.928 ± 0.034 
1.20 ±0.12 

2.100 ±0.073 
1.327 ±0.058 
0.90 ±0.12 



This table demonstrates the effects of fitting for orbital eccentricity and is to be compared with Table [6] which presents our adopted 

solutions. 

Table 6. System parameters from RV and transit data from our adopted, circular solutions 



Parameter (Unit) 



WASP-44 



WASP-45 



WASP-46 



P(d) 

T (HJD, UTC) 

T M (d) 

T12 = T 34 (d) 

a/R, 

R pl 2 /R* 2 

b 

i (°) 

K± (ms" 1 ) 
7 (ms -1 ) 

e 



rmSgpcc 

Xspcc 



(ms- 1 ) 



2.4238039 ± 0.0000087 
2 455 434.376 ± 0.00040 
0.0936 ± 0.0022 
0.0146 ±0.0026 
8.06tg;» 
0.01588 ± 0.00076 

0-56018:;™ 
86.02IJ; 11 

138.8 ± 9.0 
-4 043.82 ±0.71 
(adopted) 
25.0 
10.9 



3.1260876 ±0.0000035 
2 455 441.26925 ±0.00058 
0.0746 ± 0.0035 
> 0.0270 
9.22 ±0.74 
0.01495 ± 0.00071 
U.S»s_ 026 
84.471°;?* 

148.3 ±4.1 
4 549.95 ± 0.63 
(adopted) 
9.7 
7,1 



1.4303700 ±0.0000023 
2 455 392.31553 ± 0.00020 
0.06973 ± 0.00090 
0.0178 ±0.0013 
5.74 ±0.15 
0.02155 ± 0.00049 
0.737 ±0.019 
82.63 ±0.38 

387 ± 10 
-3 778.3 ± 1.1 
(adopted) 
34.2 
19.2 



M, (M Q ) 
Ft* (Rq) 
logs, (cgs) 

P* (pq) 
T eR (K) 
[Fe/H] 

M pl (M Jup ) 

^pl (^Jup) 

logflpl (cgs) 

Ppl (pjup) 

a (AU) 

r P i, A =o ( R ) 



0.951 ±0.034 
n O27+0-068 
u ' yz ' -0.074 



4.481_ 057 
-,+0.32 
-0.22 



i.i9: 

5410 ± 150 
0.06 ±0.10 



0.889 ±0.062 
1.14 ±0.11 

o 1Q7 +0.094 
°- xa -0.082 

61+ - 23 

u ' ul -0.15 

0.03473 ± 0.00041 
1343 ± 64 



0.909 ± 0.060 
0.945j 
4.445;! 
1.08^ 



-+0.087 
-0.071 
-+0.065 
-0.075 
i+0.27 
-0.24 



5140 ± 200 



1.007 ±0.053 
1 16+ ' 28 

1,lu -0.14 
n no + 0.11 

0.64 ±0.30 
0.04054 ± 0.00090 
1198 ±69 



0.956 ± 0.034 
0.917 ±0.028 
4.493 ± 0.023 
1.24 ±0.10 
5620 ± 160 
-0.37 ±0.13 

2.101 ±0.073 
1.310 ±0.051 
3.447 ±0.033 
0.94 ±0.11 
0.02448 ± 0.00028 
1654 ± 50 
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Figure 8. Upper panels: Periodograms for the WASP data from three seasons for WASP46. The date range (JD— 245000) is given in the 
title of each panel. Horizontal lines indicate false alarm probability levels FAP=0.1, 0.01, 0.001. Lower panels: Data folded and binned 
in 0.005 phase bins for the periods noted alongside the date range in the plot title. 



6 DISCUSSION 

WASP-44b is a 0.89- M Jup planet in a 2.42-day orbit around 
a solar metallicity ([Fe/H] = 0.06 ± 0.10) G8V star. The 
star is observable from both hemispheres and seems to be 
inactive as it shows no signs of rotational modulation or Ca 
II H+K emission. 

WASP-45b is a Jupiter-mass planet that partially tran- 
sits its metal rich ([Fe/H] = 0.36 + 0.12) K2V host star every 
3.13 days. The (near-) grazing nature of the transits led us 
to impose a prior on the density of the host star so that the 
derived stellar and planetary radii are physical. The star 
may be active as weak Ca II H+K emission is seen in its 
spectra, though no rotational modulation is evident in the 
light curves. 

WASP-46b is a 2.10-M. Tup planet in a 1.43-day orbit 
around a metal poor ([Fe/H] = -0.37 ±0.13) G6V star. We 
found the star to be active as weak Ca II H+K emission is 
present in its spectra and a clear 16-day rotational modu- 
lation signature is evident in the WASP light curves. The 
slower rotation rate suggested by spectrocopic vsinl indi- 
cates that the stellar spin axis is inclined relative to the sky 
plane (I = 41 ± 40°). By observing WASP- 46b's Rossiter- 
McLaughlin effect (e.g. ISimpson et al.l l201lh ■ which is pre- 
dicted to have an amplitude of 18 ms -1 , we could thus mea- 
sure the true angle between the stellar spin and the plane- 
tary orbit axes, rather than just the sky-projected angle. The 
large size of the planet, the high stellar insolation, and the 
relative lateness of the host star make WAS P-46b a good tar- 
get for infrared occultation photometry fe.g. lAnderson et al.l 

iboiibh . 

The upper limits placed on the lithium abundances 



suggest that each star i s at least a few Gyr old (Table [3] 
ISestito fc Randicbi r2005). which is in contrast with the gy- 
rochronological ages of ~0.9-1.4 Gyr. 

6.1 Fit for eccentricity or impose a circular orbit? 

To obtain accurate stellar densities and the derivative stel- 
lar and planetary radii, we require not only high-quality 
transit light curves, but also accurate determinations of 
the orbital eccentricities. Measurements of these parame- 
ters inform intense theoretical efforts concerning tidal cir- 
cularisation and heating, bulk planetary composition and 
the observed systematic errors in planetary and stellar radii. 
For those planets with poorly constrained eccentricity, fur- 
ther radial-velocity and occultation-timing measurements 
are necessary. However, these will be performed slowly and 
perhaps only for a fraction of systems. As such, we should 
make every effort to avoid introducing biases into the sample 
that we have to work with. 

Around half of all known transiting planets (i.e. those 
confirmed by a dynamical mass determinatio n) are hot 
Jupite rs (P < 4 d and M pJ w 0.5-2 Mj up ) and iPont et all 
|201fl ) find no evidence for an eccentric orbit among these. 
Further, tidal theory predicts that the timescale of tidal cir- 
cularisation is a sharp function of semimajor axis and planet 
radius, and a weaker func tion of planetary and stellar masses 
jGoldreich fc Soter|[l96^ . That is, low-mass, bloated plan- 
ets in short orbits around high-mass stars are expected to 
circularise quickest. So, there is a theoretical basis and an 
empirical basis for expecting a newly-discovered hot Jupiter 
to be in a circular orbit. 

When a hot Jupiter is first announced the radial- 
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velocity data often poorly constrain eccentricity and oc- 
cupation observations are almost never available. In those 
circumstances, some authors have opted to impose a cir - 
cular orbit (e.g. lAnderson et al. I l2011al . iFaedi et al.ll2011al ). 
as is theoretically and empirically reasonable, though doing 
so does artificially reduce the uncertainties on the system 
paramters . Other authors have opted to fit for an eccentric 
orbit (e.g. lHellier et al.ll20ld iBakos et al.ll201lh so as to ob- 
tain more conservative error bars. Fitting an eccentric orbit 
model to RV measurements of a (near-)cir cular orbit will 
alway s result in a non-zero value of e (e.g. lLaughlin et al.l 
2005), and this can affect the derived stellar and planetary 
parameters. As compared to the circular solution, the eccen- 
tric solution will typically change the planet's speed during 
transit and, as the transit duration is fixed by measurement, 
a change in speed necessitates a change in distance travelled 
during transit, i.e. the derived stellar radius. In turn, as the 
ratio of planetary-to-stellar radii is fixed by the measured 
transit depth, then the derived planet radius is changed. 

Being aware of these issues, some authors have opted 
to present both a circular solution, which is most likely, and 
an eccentric solution, which gi ves a reliable account of the 
uncertainties in the data (e.g. lAnderson et al.l [20101. One 
issue with this approach is that many authors tend to prefer 
to work with a single solution and so will adopt the one they 
favour regardless. Another approach is to present the best- 
fitting parameter values fr om a circular fit and the error bars 
from an eccentric fit (e.g. lEnoch et al.ll201ll . iMaxted et alj 
2011). Thus the most likely solution is presented with con- 
servative error bars, but care must be taken to calculate the 
1-a confidence intervals in the posterior probability distri- 
butions about the best-fitting circular values. 

As the best-fitting eccentricity values for the three sys- 
tems presented in this paper are small (e = 0.02-0.04) the 
differences between the best-fitting parameter values from 
circular and from eccentric solutions are small (compare 
Tables [S] and O. However, eccentricity is often poorly deter- 
mined in a discovery paper and the choice to fit for eccentric- 
ity then has a greater impact. An illustrativ e example is the 
case of WASP-I7b. In the discovery paper l|Anderson et al.l 
2010), we presented three solutions, including one with an 
imposed circular orbit which gave R p i — 1.97 ± 0.10 i?j up , 
and another, which we adopted as our preferred solution, in 
which we fit for an eccentric orbit that gave R p \ = 1-741q'23 
R 3up and e = 0.129l^. We had initially favoured the cir- 
cular solution as theory suggested this as the most likely and 
the additional two free parameters of an eccentric orbit did 
not significantly improve the fit. However, we adopted the 
eccentric solution due to its more conservative error bars. 
We had been able to place only a weak constraint on ec- 
centricity as the low mass of the planet (0.49 Mj up ) and 
the high effective surface temperature of the star (6650 K) 
caused our measurement of the reflex motion of the star 
to be quite low S/N. This situation improved when we ob- 
served two occultations of WASP-I7b, which we used to 
show the o rbit of WASP- 17b to b e very slightly eccentric 
(e « 0.0055 lAnderson et al.ll2011bh and the planet radius to 
be 1.991 ±0.081 Rj up - As such, the circular solution that we 
presented in the discovery paper was closest to describing 
the system and we would have done better to have adopted 
that. 

By the time we had shown the orbit of WASP-17b to 



be near-circular, a number of papers had, at least in part, 
based their models, results or conclusions on the values of 
one or both of e and R p \ of the discovery paper's adopted 
solut i on (e.g. S chlauf nianll 20 1 d. Hansenl2010l ICowan fe Agoll 
1201 jjEhrenreich fc DeserfeOllf T Though it is true that our 
updated solution was not excluded, not even at the 1-cr level, 
some of those studies had only considered WASP-17 due to 
its apparent significant eccentricity, and other such works 
were sure to follow. For example, the possibility that tidal 
dissipation could have inflated WASP-17b would probably 
have b een tested as was done bv llbgui. Spiegel fc Burrows! 
(|201lT ) for other systems. Before we showed the planet to be 
even more bloated than first suggested and placed a much 
tighter constraint on eccentricity, tidal heating had seemed 
a much more likely explanation for the anomalously large 
size of WASP-17b. 

A more recent example, which is similar to the case of 
WASP-17 and became available after submi ssion of this pa- 
per is that of HAT-P-32 and HAT-P-33 (|Hartman et all 
2011 !). In both systems, the signal-to-noise of the stel- 
lar radial motion, as induced by the orbiting planet, is 
low and eccentricity is poorly constrained. In an early 
draft (arXiv:1106.1212vl), eccentric solutions were pre- 
sented alongside circular solutions such that the reader could 
choose which to adopt. The derived planetary radii differ sig- 
nificantly between the eccentric and the circular solutions, 
with 1.8 Rj up versus 2.0 Rj up for HAT-P-32b and 1.7 i? Jup 
versus 1.8 Rj up for HAT-P-33b. In the published version 
of the paper, lHartman et al.l (120111 ) instead suggest adop- 
tion of the circular solutions, but caution that several ec- 
centric short-period planets are known. The three example 
exoplanets given, XO-3b, WASP-14b and HAT-P-21b, are 
each much more massive and smaller than HAT-P-32b and 
HAT-P-33b, and are in similar or longer orbits. Each effect 
acts to increase the tidal circularisation time-scale. Thus, it 
is little wonder that these planets are in eccentric orbits, but 
it would be suprising if HAT-P-32b or HAT-P-33b were. 

Unless there is clear evidence for an eccentric orbit, we 
recommend imposing a circular orbit when deriving the pa- 
rameters of a hot- Jupiter system. This will avoid creep in the 
literature of parameter values as further RV and occultation 
observations show discovery-paper eccentricities to be spu- 
rious. Though further measurements will show some planets 
to have a small, non-zero eccentricity, for the majority eccen- 
tricity will be constrained to ever-smaller values consistent 
with zero. 
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